Theory of electric polarization in multi-orbital Mott insulators 
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The interaction between the electric field, E, and spins in multi-orbital Mott insulators is studied 
theoretically. We find a generic dynamical coupling mechanism, which works for all crystal lattices 
and which does not involve relativistic effects. The general form of the coupling is — T a E a et,, where 
e is the 'internal' electric field originating from the dynamical Berry phase of electrons and T ab 
is a tensor determined by lattice symmetry. We discuss several effects of this interaction: (i) an 
unusual electron spin resonance induced by an oscillating electric field, (ii) the displacement of spin 
textures in an applied electric field, and (iii) the resonant absorption of circularly polarized light by 
Skyrmions, magnetic bubbles, and magnetic vortices. 

PACS numbers: 73.43.-f,75.70.Kw,85.75.-d,85.70.Kh,85.75.-d 
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Introduction: The manipulation of magnetic patterns 
with an applied electric field is an important issue 
both for fundamental physics and for applications to 
spintronicsi Naively, electrons in Mott insulators should 
be inert to the electric field oscillations with energies 
below the charge gap. Yet, a number of multiferroic 
materials, in which electric polarization is induced by 
spin orders^ show strong response to the electric field 
at frequencies of magnetic excitations, resulting in the 
so-called electromagnon peaks in optical absorption! 3 -! 4 - 

Two main mechanisms for the coupling between the 
electric field and spins have been identified so fax£~— One 
is the lattice and electronic polarization induced by the 
Heisenberg spin exchange energy, which gives rise to the 
'bond' electric dipoles, Py, proportional to the scalar 
products of spins: Pfj = TtfjSi ■ Sj. a The other origi- 
nates from the relativistic spin-orbit interaction inducing 
the dipole moments proportional to the vector products 
of spins, Py = ctGij x (S; x Sj), where is the unit 
vector parallel to the bond£~— The effectiveness of these 
mechanisms is restricted by symmetry requirements to 
special lattice geometries and magnetic orders, such as 
cyloidal spirals and the antiferromagnetic E-type order 
in orthorombically distorted manganites. One of the mo- 
tivations for this study is to find a generic mechanism 
that couples electric field to spin patterns in insulating 
ferromagnets independently of their crystal structure. 

Another motivation is the recent upsurge of interest 
in Mott insulators close to the transition into metallic 
state, e.g. 3d and 5d transition metal oxides,^^ organic 
crystals, 11 cold atomSfi 2 - and quantum dot arraysji 3 - The 
proximity to metallic state enhances fluctuations of the 
electron charge density playing the crucial role in the 
coupling of the low-energy spin degrees of freedom to 
the electric field. Furthermore, as we show below, large 
distances between the localized charges in artificial Mott 
insulators, 13 can strongly amplify their response to an 
applied electric field. 

In this paper, we study theoretically the electric polar- 



ization induced by time-dependent spin patterns in mag- 
netic insulators taking into account the spin dynamics 
during the exchange process. We derive the magneto- 
electric coupling for a multi-orbital Hubbard model and 
show that the most universal mechanism that does not re- 
quire special crystal lattices and relativistic effects gives 
the electric polarization proportional to the "internal" 
electric field e(x, t) associated with the Berry phase of 
dynamical spins, 



e(x,i) = (1/2) sm9(d t 6V ip - d t ifV 6). 



(1) 



This quantity has been discussed in the context of the 
"electromotive force" or "spin motive force" in metallic 
ferromagnetic systems* 1 ^— We show that this field is 
also relevant to the insulating magnets contributing to 
their dielectric response. 

The model: Our microscopic model includes the sum of 
electron Hamiltonians on transition metal sites, 



ffsite = Un(n - 1) + An/3 - J/? [S • (s Q + sp) + s a ■ sp] . 

(2) 

where the first term is the on-site Coulomb repulsion, 
n = n a +np being the total number of itinerant electrons 
on the site, A > is energy splitting between the two or- 
bitals a and j3, and the last term is Hund's rule coupling 
between the local spin Si and the spin Si a = \c\ a crci a 
of the itinerant electron on the orbital a — a, f3. The 
zeroth-order Hamiltonian is the sum of all on-site en- 
ergies and Hg describing spin interactions that do not 
originate from exchange processes. 

We do perturbation theory in the hopping energy of 
itinerant electrons, 



ia;jb 



(3) 



where tj(,,i a = t* a - 6 is the amplitude of hopping from the 
orbital a on the site i to the orbital b on the site j. In 



2 



the presence of electric field the hopping amplitudes are 
modified using the 'Peierls substitution'^ 



tjb,ia ^ tjb,ia& 



'. J Xj dx-A(x,i) 



(4) 



where A is the vector potential, — e is the electron charge 
and h = c = 1 (furthermore, the term —eA^n is added 
to the on-site Hamiltonian). 

In general, the electric field dependence of hopping 
amplitudes is more complex than the one described by 
Eq.Q. If the metal-ligand- metal bond has a nonzero 
electric dipole moment, an applied electric field paral- 
lel to the dipole results in a linear dependence of the 
hopping amplitude on E, which gives rise to the ex- 
change striction and relativistic mechanisms discussed in 
the introduction^ - — Here we only take into account the 
electric field dependence resulting from the Peierls sub- 
stitution present for any bond geometry. 

The model with two orbitals per site is the simplest 
model of a multi-orbital Mott insulator. We first show 
that the electron hopping between different orbitals, e.g. 
the orbital a on site 1 and the orbital j3 on site 2 
(see Fig.QJ, favoring ferromagnetic interactions between 
spins, gives rise to a dynamical electric polarization in- 
duced by rotating spins. Then we explain why this does 
not happen in the single-orbital Hubbard model. For sim- 
plicity we assume that there is only one itinerant electron 
per site and consider the strong Hund's rule coupling and 
large 5* limit, in which the spin of the itinerant electron 
is parallel to the local spin. 

Electric polarization of spins: First we discuss interac- 
tions between two sites, 1 and 2, in the ferromagnetic 



model, in which only tj a 



t 



t ^ (see Fig. P. 



To second order in the hopping amplitude t the correction 
to the imaginary time spin action is 



5S = t 2 J dn J dT ie - u to" 7 *) [Caicai + C 12 c 12 } , 



(5) 



p — — p 



(a) 




FIG. 1: (a) Two-orbital model of electrons interacting with 
the local spins Si and Sj through the Hund's rule coupling Jh- 
(b) Physical realization of the two-orbital model describing 
the hopping between the a = d yz and ft = d xz orbitals of the 
magnetic ions on the metal-oxygen plaquette mediated by the 
p z -orbitals of oxygen ions. 



where U' = U + A- J H /A, 

-ie f ! dr J" 2 dx-E 



C21 



C21 



(6) 



P dr[<n 1 (r)|a T |n 1 (r)>~<n 2 (r)|a T |n 2 (r))-iH 21 ] 



x(ni(Tf)|n 2 (Tf))(n 2 (Ti)|ni(Ti)), 



with 5H21 
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2(s+i/2) \^ ' aSI - ni ' SSI 
the spin energy change due the electron hopping from 
site 1 to site 2 (C12 , c\ 2 are obtained by interchanging 
the indices I and 2). Here, n^a is the unit vector in 
the direction of the local spin Si^ and |n) denotes the 
eigenstate of the itinerant electron with spin parallel to 
n. 

C21 is invariant under local gauge transformations and 
so is C21 • Introducing the vector potential of the 'internal' 
field by ao = i(n\d T \n) and a = i(n|<9 x |n), we can write 
the overlap of the electron spin wave functions in the 
form, 



Hs describing 



njni 



"21 

cos — e 
2 



1 eZx-a 



(7) 



where #21 is the angle between ni and n 2 . The vector 
n(x) that defines a varies between ni to n 2 along the 
segment of a circle on the unit sphere, while x varies 
from xi to X2. We now can write C21 in the manifestly 
covariant form 



C21 



921 (Tf) 02l(Ti) pdT(ip 

— - cos -^-e J ^ Jx i 

2 2 



dx-G — SHi2 



(8) 



where e = <3 x ao ~ d T a. is the gauge invariant internal 
electric field. Comparing Eqs. ([()]) and ©, we find that 
the correction to action only depends on the combination 
of applied and internal electric fields, eE — e. 

Since the time spent by the hopping electron on a 
neighboring site, Tf — n ~ (U')^ 1 , is much shorter than 
the characteristic time of spin dynamics, Cy and in 
Eq.© can be expanded in powers of 71 — 71, which gener- 
ates an expansion of the spin action in powers of (U')~ l . 
To lowest order we obtain an effective ferromagnetic in- 
teraction between the spins, 



v 



'ni ■ n 2 + 1) 



(9) 



The third-order term in the expansion gives the interac- 
tion described by the real time Lagrangian, 



L = 



t 2 



(U>) 



(n r n 2 + 1) 



fX 2 

/ dx-(eE-e) + SH 21 

Vxi 



(10) 

The term oc E 2 is the spin contribution to the dielectric 
susceptibility, while the term linear in E describes the 
coupling of the external electric field to the spin-induced 
electric polarization: 

L E = £ {(hi + ri 2 ) • [ni x n 2 ] + 2 (ni • n 2 + 1) SH 2 i} , 

(11) 
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where £ 



feE-(x2-xi) 
(U< 



The second term in curly brack- 



ets describes the charge re-distribution between the two 
sites that takes place when ni ^ 112 , while the first term 
is the dynamical polarization originating from the cou- 
pling between the external and internal electric fields, 
E ■ e, which for insulators was not considered before. 
Though weak, it leads to a number of unusual effects 
discussed below. 

Physical consequences: The dynamical part of Le is elim- 
inated by the rotation of the spins around [ni x 112] , 



5xii 
Sn 2 



T72) ( n 2 - ( n i ' n 2)ni) , 



(S+l/2) 



" (s+1/2) ( n i ~ (ni-n 2 )n 2 ), 



(12) 



applied to the Berry-phase term in the spin Lagrangian, 
£s = 0S+i) ^ (cos 0,-1)^, (13) 



i=l,2 



where Bi and fi are the polar angles describing n^. If the 
spin Hamiltonian, Hg, is not rotationally invariant, this 
transformation generates a coupling to electric field in the 
spin Hamiltonian. For example, the magnetic anisotropy, 



Ha 



A 
~2 



(sif + (§£)' 



where Si = Si + Si a , gives rise to the interaction 



-Hint — 



A£(t) 
S 



(14) 



(15) 



(S = S + \). The time-dependent electric field induces 
transitions between the states of opposite parity under 
the permutation of the two sites, corresponding to even 
and odd values of the total spin of the two sites, St, 
while the projection of the total spin on the anisotropy 
axis, Sf, remains constant. These selection rules have to 
be compared with those for the conventional ESR where 
the time-dependent magnetic field leaves the total spin 
unchanged, while its projection on the static magnetic 
field changes by ±1. 

Consider now an arbitrary spin texture in a ferromag- 
net with n(x, t) varying slowly at the lattice constant 
scale. The Lagrangian describing the linear coupling of 
the texture to the electric field is given by 



L E = - Jd 3 xT ab E a e b (x,t), where 



(16) 



T 



ah 



e 1 



{U'f v 



(v is the unit cell volume). For a simple cubic lattice 
with the lattice constant a, T ab = -^5 ab with g = ^jr- 
Importantly, Eq. (|16[) entirely comes from the dynamical 



part of Eq. (ITT|) . as the static polarization in the contin- 
uum limit is a total derivative. 

Similarly to the two-spin case, the transformation 
n(x, t) ->■ n(x + X,i) with X = ^E, applied to the 
Berry phase term, cancels the interaction Eg. ([TBI). Since 
this transformation leaves the Hamiltonian of a transla- 
tionally invariant system unchanged, the effect of electric 
field is to shift a spin texture as a whole by the vector 
—X. The shift is a small fraction of the lattice constant: 
"a" ~ \u^)* ' ^ can ' h° wever i be enhanced by proxim- 
ity to the metal-insulator transition (through larger t/U 
ratio) and by a larger distance between the spins, a, in 
magnetic semiconductors and quantum dot arrays. 

A much stronger effect of the interaction Eq. flTB")) is 
the resonant absorption of circularly polarized light by 
Skyrmions, magnetic bubbles, and magnetic vortices. 
These magnetic defects in two spatial dimensions carry a 
nonzero topological charge^ Q = j- Jd 2 xn ■ d x ii x d y n, 
integer for Skyrmions/bubbles and half- integer for vor- 
tices. 

Magnetic vortices are spontaneously induced by mag- 
netostatic interactions in nanodiscs of ferromagnetic 
metals, 19 while periodic arrays of magnetic bubbles ap- 
pear in thin-film ferromagnets with a strong out-of-plane 
anisotropy upon application of magnetic field on the or- 
der of 10 2 Oej 2 ^ Similar arrays of skyrmions, which are 
bubbles with "thick" domain walls, were recently ob- 
served in bulk ferromagnetic metals without inversion 
symmetry i 22 i 23 

According to Eq. fToT) . a moving topological defect in- 
duces the net electric dipole moment in the direction nor- 



mal to its velocity, P oc gQ 



x R 



where R = (R x , R y ) 



is the position of the center of the defect and z is the 
unit vector normal to the film. Consider such a defect 
in a ferromagnetic insulator with a confined geometry, 
which breaks translational symmetry, e.g. a vortex in a 
nanodisc. We assume that the confining potential has 
the form, U = ^{R? x + Ry). In the adiabatic limit the 
dynamics of the collective coordinates Ri is described by 
Thicle equations^— 



G, 



Rj + ^Ej 
3 s j 



aTijRj 



dU 
dR~ 



(17) 



to which we added the coupling of spins to the electric 
field. Here, a is the Gilbert damping constant and the 
nonzero components of the tensors G%j and Tij are, 



Gxy — —G yx — 4nQ and 
Fxx = T yy = Jd 2 xdin ■ dill. 



(18) 



In absence of electric field the (damped) eigenmode 
R(t) cx (cos57t, — q sin fit) describes the rotational mo- 
tion of the center of the spin texture with the frequency 
fi = in the direction defined by q = sign(Q) = ±1. 

The response to the rotating electric field 



E(i) = E u (cosujt, —crsinwt) . 



(19) 
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(cr = ±1) at the resonant frequency, u — f2, is given by 

Xn = 9 -§^{ (20) 
b I -2^H7' for a =- g , 

where r = is the relaxation time. For fir ~ ^ < 1 
the excitation of the rotational motion by the electric 
field with a = +q is resonantly enhanced by the fac- 
tor ^ compared to the shift in translationally invari- 
ant systems, while for a = —q there is no enhancement. 
For magnetic insulators with a ~ 10~ 3 — 10~ 2 (see e.g. 
Refl20|) the resonant enhancement by three to four or- 
ders of magnitude opens a possibility to manipulate spin 
textures with an electric field. 

Discussion: Now we discuss why the presence of several 
orbitals is essential to obtain the linear coupling of spins 
to electric field Eq. (fTTj) and why such a coupling does not 
exist in the single-orbital model, in which only tj a ^ a ^ 
and interactions between spins are antifcrromagnetic. In 
that case 

P dr[(n 1 (r)|a T |n 1 (r)>-<-n 2 (r)|a T |-n 2 (r))] 

c 2 i = e J n 

x( ni (Tf)| -n 2 (7t)>(-n 2 (7i)|ni(7i)) (21) 

(we ignore the changes in the magnetic energy SH21, 
since the static polarization disappears in the contin- 
uum limit). We note that for a particular choice of 
spinor wave functions, describing electrons with spin par- 
allel/antiparallel to the local spin, 

'-> = (~ttf )■ « 

C21 is symmetric with respect to the interchange of the 
indices 1 and 2: en = ci\. Since the E-dependent term, 
C21, is antisymmetric with respect to this permutation, 
there is no linear coupling of spins to electric field. Due 
to the gauge invariance of C21 this result is independent 
of choice of spinor wave functions. 

This result can be physically explained as follows. The 
first term in Eq. (|ll[) originates from the spin dynamics in 
the virtual states with two electrons occupying the same 



site. More precisely, it describes the difference between 
the spin rotations in the state where both electrons oc- 
cupy site 1 and the state with two electrons on site 2. In 
the single-orbital case, however, the virtual state of two 
itinerant electrons is a spin-singlet independent of which 
site is doubly occupied and which is empty. Hence, no 
dynamical linear coupling to electric field. 

In multi-orbital Mott insulators both exchange pro- 
cesses described by the models I and II take place. Im- 
portantly, the coupling of spins to the electric field, ex- 
pressed by Eqs. and (|16|) . is present even if the dom- 
inant interaction between spins is antiferromagnetic, i.e. 
our results apply to insulators with ferrimagnetic, canted 
and spiral spin orders. 

There are several interesting problems left for future 
studies. One is the relevance of the present mechanism 
for spin liquid states,— where effects of the spin Berry 
phase are enhanced by strong spin fluctuations. Another 
issue of interest is the ring-exchange processes giving rise 
to persistent orbital currents in Mott insulators^ Dy- 
namical effects resulting from the ring-exchange process 
deserve a scrutiny. 

In summary, we show that the dynamical spin Berry 
phase in multi-orbital Mott insulators couples the elec- 
tric field to the translational modes of spin textures. We 
derive equations of motion for the center-of-mass coordi- 
nates of Skyrmions and magnetic vortices in an applied 
electric field and predict the resonant absorption of the 
circularly polarized light by these topological objects as 
well as the ESR effect where spin transitions are induced 
by the time-dependent electric field. 
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